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We study dynamical recurrences of a Bose- Einstein condensate in an optical crystal subject to a 
periodic external driving force. The recurrence behavior of the condensate is analyzed as a func- 
tion of time close to nonlinear resonances occurring in the classical counterpart. Our mathematical 
formalism for the recurrence time scales is presented as delicate recurrences, which take place, for 
instance when the lattice is perturbed weakly and robust recurrences, which may manifest them- 
selves for a sufficiently strong external driving force. The analysis not only is valid for a dilute 
condensate, but also is applicable for strongly interacting homogeneous condensate provided; the 
external modulation causes no significant change in density profile of the condensate. We explain 
parametric dependence of the dynamical recurrence times which can easily be realized in laboratory 
experiments. In addition, we find a good agreement between the obtained analytical results and 
numerical calculations. 

PACS numbers: 67.85.Hj, 42.65.Sf, 03.75.-b 



I. INTRODUCTION 

Coherent control of matter waves is at the heart of 
many theoretical and experimental interests. Coherence 
of matter waves in optical crystals has evolved as an ac- 
tive research area over last two decades and led to an- 
alyze fundamental characteristics of quantum mechanics 
[1], quantum tunneling [2], Bloch oscillations [3], quan- 
tum information [4] and quantum chaos[5, 6]. The fas- 
cinating developments in this subject have wide applica- 
tions from quantum metrology [7] to quantum corrals [8] . 
Time periodic modulations in matter wave optics have 
given birth both to hybrid nano-opto-mechanical sys- 
tems [9] and driven billiards [10]. Hence, the theoretical 
and experimental advancements have led to explore co- 
herent transport [11, 12], controlling ratchet effects [13], 
dynamical localization [14, 15], coherent destruction of 
tunneling [16], photon assisted tunneling [17], entangle- 
ment [18], precision measurement of gravitational accel- 
eration [19] and extracting nearest-neighbor spin correla- 
tions in a fermionic Mott insulator [20] in driven optical 
crystals. Furthermore, realization of strongly correlated 
phases [21], coherent acceleration of matter wave [22] in 
driven optical crystals [23] have been studied in great 
details. 

A quantum particle in its early evolution in a bounded 
system follows classical mechanics and reappears after a 
classical period following classical trajectory. Later, fol- 
lowing wave mechanics, it spreads and collapses as a con- 
sequence of non-linearity of the energy spectrum. How- 
ever, the discreteness of the quantum mechanics leads 
to its reconstruction at various longer time scales, which 



define quantum recurrence times, such as, quantum re- 
vival time and supper revival time. We explain these 
recurrences of a quantum particle in a quantum chaotic 
system, which manifests complex dynamics in its classi- 
cal counterpart. Due to their dependence on modulation 
effects, these dynamical recurrences are different from 
those taking place in an undrivcn system, and may be 
considered as a probe to study quantum chaos [24]. Dy- 
namical recurrences, in a bounded driven system, origi- 
nate from the simultaneous excitation of discrete quasi- 
energy states [25, 26]. In the present paper we explain 
matter waves dynamics in a phase modulated optical 
crystal and provide analytical relations for classical pe- 
riod, quantum revival and super-revival time scales. In 
addition, we explain the existence of delicate dynamical 
recurrences, for weak effective modulation, and robust 
dynamical recurrences for strong effective modulation 
of the periodically driven optical crystal. The classical 
counterpart of the dynamical system displays dominant 
regular dynamics and dominant stochastic dynamics, one 
after the other, as a function of increasing modulation 
amplitude [6]. Our analysis is valid for weakly interact- 
ing BECs, a situation that can be realized experimentally 
with Feshbach scattering resonances [27] . In addition, we 
suggest that the analysis is valid for strongly interacting 
homogeneous condensates as well, where, nonlinear term 
can be replaced by an effective potential provided the ex- 
ternal modulation causes slight changes in density profile 
of the condensate as discussed in Ref. [28]. Hence, due 
to spatial and temporal periodicity in the driven optical 
crystal, the quantum dynamics of the condensate inside a 
nonlinear resonance is effectively mapped on the Math- 
ieu equation. Our analytical findings are confirmed by 
numerical results. 
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The paper is organized as follows: In Sec. II, we 
model the effective Hamiltonian governing the time evo- 
lution. Later, the quasi-energy spectrum for nonlinear 
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resonances in periodically driven optical crystal is ex- 
plained. We obtain mathematical relations for recur- 
rence time scales for the two cases; delicate dynamical 
recurrences and robust dynamical recurrences for matter 
waves in driven optical crystal in Sec. III. The analytical 
and numerical results arc summarized in Sec. IV. 



For the dynamics of weakly interacting BECs in optical 
crystal the interaction term can safely be neglected [14], 
which is a situation achievable in the present day experi- 
ments with the advent of Feshbach scattering resonances. 
In addition, the effective potential ^ cos 2z -t- G|-0p seen 
by each atom may as well be written [28] as, 



II. THE MODEL 



V 

Veff ~ — cos(2z) -I- const, 



We consider a condensate which is loaded in the lower 
edge of the band corresponding to optical crystal to avoid 
dynamical and Landau instabilities [29]. The dynamics 
of the matter wave in driven one dimensional optical crys- 
tal, strongly confined by radial trap is governed by the 
Hamiltonian [15, 30], 



ALsin(w„t)}] + .giD|V'|^, 



1+4G- 



where, V 

The analytical result, calculated using perturbation 
theory [28], is valid as long as the condensate density 
is nearly uniform, i.e. , << 1 which describes either 
a weak external potential VJ, or a strong atomic interac- 
tion G. The condition is experimentally confirmed for 
one dimensional potential [31]. The effective dynamics 
of matter wave in this regime is governed by the Hamil- 
tonian 



where, is wave number, Vo defines the potential depth 
of the lattice. Moreover, AL and ujm ai'c respectively, am- 
plitude and frequency of external drive, whereas, M is the 
mass of an atom. Furthermore, guj = hki^u^_ag defines 
the effective two body interaction coefficient, uj^_ is radial 
trap frequency and is inter-atomic s-wave scattering 
length. 

The unitary transformation^ ?/; = 
'!/;(a-,t)cxp[^{a;„iMALcos(wmt)a; -|~ /3(t)}], (where, 

cj^jAL^A/ rsin(2w„t) 



m 



t\), for a frame co- moving 



with the lattice, modifies the Hamiltonian as 



2M 2 



FxsiiiLL>mt + gmlM'^- (1) 



Here, F = M ALw^ is amplitude of inertial force emerg- 
ing in the oscillating frame. To examine the dynamics of 
cold atoms in driven optical crystals numerically, Hamil- 
tonian (1) is expressed in dimensionless quantities. We 
scale the quantities so that z = kj^x, r = oJmt, "0 = 
where, ^/r^o is average density of a condensate. Mul- 



tiplying the Schrodingcr wave equation by where 



2M 



is single photon recoil frequency, we get dimen- 



sionless Hamiltonian, viz 

k^ 92 



H 



2 dz 



H cos2z + Azsinr + GJV'I^ 



(2) 



Here. 



and A = 



F cIlT: 

huJm 



= ki^AL is scaled 



modulation amplitude. Also Vo 



2hu]„ 



T is scaled time 



in the units of modulation frequency, Wm, and de- 



fines the periodicity of optical crystal, 
rescaled Planck's constant is = — . 



In this case, the 



^2 q2 ^ 

H = - — —^ + — cos2z + XzsinT. (3) 
2 oz^ I 

Here and onward V is expressed in the units of recoil 
energy i?,, hw^. 

In order to understand the classical dynamics of the 
system, we plot the Poincare surface of section for shal- 
low optical crystal iy = 2) with modulation strengths 
A = 0, 0.5, 1 and for deep lattice {V ~ 16) with mod- 
ulation strengths A = 0, 0.5, 1.5 as shown in Fig. 1. 
From phase space plot, we see the appearance of 1:1 res- 
onance for A > 0. This resonance emerges when the time 
period of external force matches with period of unper- 
turbed system. One effect of an external modulation is 
the development of stochastic region near the separatrix. 
As modulation is increased, while the frequency is fixed, 
the size of stochastic region increases at the cost of reg- 
ular region. 




^ The unitary transformation is time periodic and preserves the 
quasi-cncrgy spectrum. 



FIG. 1: Poincare phase space for different modulation 
strengths and lattice potentials. Upper row: A = 0, (a) , 0.5, 
(b) 1.5, (c) and V = 2\ Lower row: A = 0, (d) 0.5, (e) 1.5, (f) 
and V = 16. 
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In order to solve the time dependent Schrodinger equa- 
tion near the nonlinear resonances of the driven quantum 
system corresponding to the Hamiltonian (2), we [32-34], 
provide an ansatz as 



where, En is the mean energy, C„(t) is time dependent 
probability amplitude, n is mean quantum number, |n) 
are eigen states of undriven system and N is resonance 
number. 

The Floquet energy spectrum for nonlinear resonances 
is given as [35, 36], 



f 
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mod Tiu!„ 



(5) 



where, a^{ii{j), q) is Mathieu characteristics parameter 
and g = is the effective modulation. Here, /3o = , 
and V =< n|z|n±l > are matrix elements and are 
considered constant at potential minima. The solu- 
tions are defined by 7r-periodic Floquet states, i. e. 
Ptj.{(^) = -f^C^ + ^^'^ A* is the characteristic exponent. 
In order to have 7r-periodic solutions, we require /i to be 
defined as /i = = 2j/N, where, j = 0, 1, 2, ....,N— 1. 

The allowed values of fj.{j) can exist as a character- 
istic exponent of solution to the Mathieu equation for 
discrete u (which takes integer values) only for certain 
value a,^{p{j),q), when q is fixed. The index v takes the 
definition v = ^^"^^"^ . Moreover, a defines the winding 
number. For driven optical potential, keeping the peri- 
odicity of Floquet solutions [34] , we take only even values 
of the index [37]. Hence, the re-scaled Mathieu charac- 
teristic exponent becomes 



iy = 2{l + 13), where, /3 



(6) 



Comparing the coefficients of eigen energy of the un- 
driven system and equation of motion of probability am- 
plitude in the absence of modulation, we get I = 
which is new band index for nonlinear resonance and at 
the center of resonance 1 = 0. Analytical expressions of 
a^{^{j),q) for two asymptotic cases i.e. for g ^ 1 and 
g ^ 1 are different [37]. We discuss these two cases in 
Sec-III separately. 



governed by the Hamiltonian 



H = 



V 



■ cos 22 



(7) 



The corresponding time independent Schrodinger equa- 
tion is eventually a Mathieu equation, where, go = -jj 
is effective crystal potential scaled by recoil energy. In 
shallow lattice potential limit, i.e. qo ^ 1, the classi- 



cal frequency, w 
C = 2 



2fi{l 



2{f. 



r} and non-linearity. 



(s2_"^)3 ■ The classical period of the conden- 
sate is inversely proportional to the average classical fre- 
quency, cj, and calculated [38] as 



id) 



2(n2 - 1)2 



(8) 



In the long time evolution the condensate spreads and 
observes a collapse as a result of destructive interference 
of constituent wavelets. At longer time the construc- 
tive interference overcomes and leads to reconstruction 
of the condensate at quantum revival time, that is in- 
versely proportional to the nonlinearity of the spectrum 
present around the mean quantum number, n, viz. 



{rev) 



27r <^ 1 



1) 



2 (n2 - 1)3 

Whereas, the super-revival time, obtained as, 
(spr) _ 7r(n2 - 1)4 



qln{v?- + 1) ' 



(9) 



(10) 



is inversely proportional to incremental change in non- 
linearity with respect to quantum number n, at mean 
quantum number n. 

In deep optical potential limit, i.e. qo » 1, the clas- 

g -) and non-linearity 

In this case, the classical time 



sical frequency is = 4(y^ — 2" 

^ _ 3(2ra+l 



21 



is C = 

period, quantum revival time and super revival time are, 
respectively, calculated as 



id) 



2v^ 



1 



/qo 



3(^2 + 1) 
28go 



(rev) 



47r 1 - 



3s 
16qo 



(11) 



(12) 



III. RECURRENCE TIMES IN A DRIVEN 
OPTICAL LATTICE 

As recurrence time scales of driven optical crystal are 
expressed in terms of undriven system, we present classi- 
cal period, quantum revival time and super revival time 
scales of matter wave in undriven optical crystal [38, 39]. 

In the absence of external forcing, the crystal poten- 
tial replaced with effective potential in equation (1) is 



and T(5'^''' = 327r^, (13) 

where, s = 2n -I- 1. 

In the following we explain the quantum recurrences of 
condensate in optical crystal in the presence of external 
periodic forcing. For the sake of clarity and simplicity, 
we classify our later discussion as delicate dynamical re- 
currences for g ^ 1 and robust dynamical recurrences for 
q » 1. 
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FIG. 2: Left panel: Classical time period versus A = ^ for 
weak external modulation (a) and for strong external modula- 

tion (b), where, /3o = for = 1. Middle panel: Quantum 
revival time versus A for weak external modulation (c) and for 
strong external modulation (d). Right panel: Super revival 
time versus A when external modulation is weak (e) and for 
strong external modulation (f). In this figure, for deep lattice 
V — 16, for shallow lattice V — 2 and 'k — 0.5. 



A. Delicate dynamical recurrences 



For weakly driven, shallow or deep, lattice (g ^ 1) 



the recurrence time scales, T^'''^ T^''"'"'' and T^"^'' 
primary resonance iV = 1 are calculated as [35] 



(cl) 



{cl) 



2 {4(/ + /?)2-l}2 



A, 



for 



(14) 



rp{rev) _ rp(rev) 



12(/ + /3)2 



2 {4(Z + /3)2-l}3 



and T 



where, A 



(spr) 
A 



2C*92(/ + /3){4(? + /3)2 + i}' 



(15) 



(16) 



(1 



■^)^^ ■ Here, time scales for weakly 
driven shallow lattice and weakly driven deep lattice are 
similar in structure, however, they differ due to different 
energy spectra corresponding to undriven crystal in the 
two cases which leads to different undriven time scales. 

In order to observe the dynamics of a condensate inside 
a resonance, we evolve a well localized Gaussian conden- 
sate in the driven optical crystal. Numerical results are 
obtained by placing a condensate in a primary resonance 
with iV = 1. 

The classical period time, quantum revival time and 
super revival time of matter waves in modulated optical 
crystal near nonlinear resonances versus scaled modula- 
tion A = ^ is shown in Fig. 2. In each plot of the figure. 




FIG. 3: Spatio-temporal behavior of atomic condensate for 
g = 0(A = 0), 1/ = 2 (a) q = 0.2^o, V = 2(b) and q = 0.2/3o, 
V — 0.3 (c). Other parameters are Ap = 0.5, /5 
S = 1. 



TV and 



left vertical axis shows the time scale when shallow op- 
tical crystal is modulated and right axis shows the time 
scale when deep optical crystal is modulated. The up- 
per row of Fig. 2 represents the time scales for small q 
values i.e. delicate dynamical recurrences as a function 
of modulation A, whereas, lower row represents the time 
scales when (7 3> 1, as a function of modulation A, i.e. 
robust dynamical recurrences. Here, left column shows 
the results related to the classical periods. Quantum re- 
vival times are plotted in middle column, whereas, right 
column shows super revival times. 

We note that when optical crystal is perturbed by weak 
periodic force, the classical period increases with modu- 
lation, as given in equation (14). Classical period for 
weakly driven shallow lattice potential changes slowly 
as compared to weakly driven deep lattice potential as 
shown in Fig. 2a. Quantum revival time in nonlinear res- 
onances versus modulation is shown in middle column of 
Fig. 2. For delicate dynamical recurrences (Fig. 2c), the 
quantum revival time decreases as modulation increases. 
The behavior of quantum revival time is given by equa- 
tion (15) for delicate dynamical recurrence. For weakly 
driven shallow optical crystal or weakly driven deep lat- 
tice, qualitative and quantitative behavior of revival time 
is almost similar. Classical period and quantum revival 
time for delicate dynamical recurrences show good nu- 
merical and analytical resemblance for the system with 
our previous work [25, 26, 40]. 

Fig. 3 shows spatio-temporal evolution of an initially 
well localized condensate in a crystal potential well. Fig. 
3a is spatio-temporal dynamics in undriven crystal, while 
Fig. 3b & Fig. 3c present the case, for external modula- 
tion q = 0.2/3o and different value of V. Spatio-temporal 
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evolution of the condensate in optical crystal shows that 
condensate diffuses to the neighboring lattice sites by 
tunneling and splits into smaller wavelets. Later, these 
wavelets constructively interfere and condensate revival 
takes place. Revival time calculated numerically is the 
same as obtained from analytical results. Keeping mod- 
ulation constant as g = 0.2/3o but for different V, which 
may be a consequence of varied atom-atom interaction, 
revival time is modified. We note that in the absence of 
interaction term, G, revival time changes with V and in- 
terference pattern is similar. But with the introduction 
of interaction term not only revival time is modified due 
to change in V, interference pattern is also modified as 
seen in Fig. 3c. 



B. Robust dynamical recurrences 

On the other hand, for strongly driven optical crystal, 
the effective modulation, q ^ 1. Time scales for the 
atomic condensate for primary resonance with = 1 
are given as [35] 



T. 



id) 



2n 



(rev) _ »7r 



3{4(/ + /3) + l} 

i6Vg 



(17) 



(18) 



crystal is given by equation (17). The behavior of classi- 
cal period for strongly driven lattice versus modulation is 
qualitatively of the same order for both strongly driven 
shallow lattice and strongly driven deep lattice as shown 
in Fig. 2b. The behavior of classical period in strongly 
driven crystal case is understandable as strong modula- 
tion influence more energy bands of undriven crystal to 
follow the external frequency and near the center of non- 
linear resonance the energy spectrum is almost linear [35] 
with assumptions q 3> 1 and / is small, i.e. condensate is 
placed near the center of resonance. 
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FIG. 4: Spatiotemporal behavior of atomic condensate. Wave 
packet dynamics in the absence of modulation (a), and in the 
presence of modulation when q = 3/3o (b). Other parameters 
are V = 0.36, Ap = 0.1 and S = 0.16 



and 



-,{spr) _ 327r 



(19) 



In case of strongly driven crystal, when external mod- 
ulation frequency is close to the harmonic frequency, ma- 
trix elements, V can be approximated by those of har- 
monic oscillator and effective modulation, q, can be ap- 



proximated as q 
tion the time scales are 



Under this approxima- 



T. 



167r(7l/VC 



16(n + 1)3 Va - {4(/ + 13) + IjqSkVC 



(20) 



t: 



(rev) 
X 



Stt 



1 - 



3{4{l + (3) + l}qSk^ 



and 



t: 



(spr) 



32{n + l)iVX 



647r(n + 1)3\/A 

3 I ■ 

k'^(2q§ 



(21) 



When lattice is strongly modulated by an external pe- 
riodic force, the classical period decreases as modulation 
increases. Classical period for strongly driven optical 



Quantum revival time in nonlinear resonances versus 
modulation is shown in middle column of Fig. 2. For 
robust dynamical recurrences, the behavior of quantum 
revival time is given by equation (18). The qualitative 
behavior of revival time for strongly driven shallow lat- 
tice is different from that of strongly driven deep lattice 
Fig. 2d, as in the later case, change in revival time is al- 
most one order of magnitude larger than the former case, 
for equal changes in modulation. Here, the difference in 
the revival time behavior for strongly driven shallow lat- 
tice as compared to strongly driven deep lattice is due 
to the difference in energy spectrum of undriven system. 
In deep lattice, due to small non-linearity more energy 
bands are influenced by the external drive and resonance 
spectrum is similar to that of harmonic oscillator near 
the center of nonlinear resonance. As modulation is in- 
creased more and more energy bands are influenced by 
external drive. Fig. 4 shows spatio-temporal evolution of 
an initially well localized wave packet in a lattice poten- 
tial well inside a resonance for F = 16. Fig. 4a is for the 
spatio-temporal dynamics of atomic condensate in the 
absence of periodic modulation, while Fig. 4b presents 
the case when external modulation, q = 3/3o- The quan- 
tum revival times in Fig. 4 seen numerically are the same 
as obtained analytically in equation (18). 
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The behavior of super revival time versus modulation 
is shown in right column of Fig. 2. For robust dynamical 
recurrences, equation (19) gives the time scale for super 
revivals. The super revival time for robust dynamical 
recurrences increases with modulation, as shown in the 
Fig. 2f. Here, qualitative behavior of super revival time 
is same for strongly driven shallow lattice and strongly 
driven deep lattice but quantitatively super revival time 
increases almost two times faster. 

Square of auto-correlation function (| A{t) |^) for the 
minimum uncertainty condensate is plotted as a function 
of time in Fig. 5 and Fig. 6 for g = 0.5/?o and q = 
1.5/3o respectively. Other parameters are Ap = Az = 0.5, 
k =0.5 and V ~ 16. In these figures upper inset is 
an enlarged view which displays classical periods, while 
lower inset of the figures show quantum revivals. Lower 
panel shows the existence of super revivals. The classical 
period, quantum revival time and super revival time arc 
indicated by arrows and showing there characteristics. 
Numerical results are in good agreement with analytical 
expressions. 




FIG. 5: Square of auto-correlation function of Gaussian con- 
densate plotted as a function of scaled evolution time. Here, 
q = 0.5/3o, S = 0.5, V = 16, Az = Ap — 0.5 and condensate 
is initially, well localized at the central lattice well around the 
lowest band of undriven crystal. Classical period, revival time 
and super revival time are indicated by arrows. 



IV. DISCUSSION 

We presented above a discussion on the occurrence of 
delicate and robust dynamical recurrences of matter wave 
in optical crystal in the presence of periodic modulation. 
The delicate dynamics condition, i. e., g < 1 is satisfied 
when shallow or deep potential is slightly perturbed by 
small external periodic force. Classical periods of a con- 
densate initially localized near the center of resonance 




FIG. 6: Square of auto-correlation function of Gaussian con- 
densate is plotted versus scaled time. Here, q — 1.5/3o, other 
parameters and conditions are same as in Fig. 5 

increases with modulation for delicate dynamics. While, 
quantum revivals and super revivals decrease with modu- 
lation. Similarly, condition g ^ 1 can be satisfied for ex- 
ample when shallow or deep potential is strongly modu- 
lated. Classical periods in robust case decrease and super 
revival times increase with increasing modulation. Here, 
quantum revival time for the case when deep lattice is 
strongly modulated increases with modulation, while, it 
decreases when shallow is strongly modulated. The dif- 
ference in the behavior is due to the contrast in energy 
spectrum of undriven crystals. When modulation is in- 
creased in deep lattice case more and more energy lev- 
els are influenced by external modulation (Fig. Ic, If), 
non-linearity near the center of resonance decreases and 
revival time increases. 

Parametric dependence of analytical results are con- 
firmed by exact numerical solutions both for delicate 
and robust dynamical recurrences. Both spatio-temporal 
(Fig. 3 and Fig. 4) and temporal (Fig. 5 and Fig. 6) 
dynamics show that the non-linearity of the un-driven 
system, and the initial conditions on the excitation con- 
tribute to the quantum revival time and super revival 
time. 

The suggested theoretical results may be realized in 
experimental set up of recently performed experiments 
at Pisa [41], where, dynamical control of matter wave 
tunneling is studied in strongly shaken crystal. A BECs 
of about 5 X 10"' ®^Rb atoms was evolved in a dipole trap 
which was realized using two intersecting Gaussian laser 
beams at 1030 nm wavelength and a power of around 1 W 
per beam focused to waists of 50 fim. After obtaining 
pure condensates trap beams were readjusted to obtain 
elongated condensates with the trap frequencies (80 Hz in 
radial and « 20 Hz in the longitudinal direction). Along 
the axis of one of the dipole trap beams a one-dimensional 
optical crystal potential was introduced and the power of 
the lattice beams ramped up in 50 ms in order to avoid 
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excitations of the BEC. The optical crystals used was cre- 
ated using two counter-propagating Gaussian laser beams 
{Xl = 852 nm) with 120 fj,m waist and a resulting optical 
crystal spacing = Al/2 = 0.426 /im. The depth Vo 
of the resulting periodic potential is measured in units 
of Er = h'^ / {2M d\) . In laboratory, accessible scaled 
optical lattice depth V ranges from 1 to 20. 

For optical crystal with potential depth V — 2, and 
V = \Q, the mean separation of the two lowest bands is 
« i.lbkHz and 20.784fci7z, respectively. For the driv- 



ing frequency, a;m/27r, ranging from 5kHz — 9kHz, the 
rescaled Planck's constant k ranges from 0.668 to 2.066. 
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